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£>. . Abstract 

(N ' 

We prove two stability results for the scale invariant solutions of the nonlinear 
heat equation dtu = Au— \u\ p ~ l u with 1 < p < 1+^, n being the spatial dimension. 
The first result is that a small perturbation of a scale invariant solution vanishes 
as t — > co. The second result is global, with a positivity condition on the initial 
■ data. 



1 Introduction 

o 

f3 \ In Q, we investigated the long-time asymptotics of the solutions of nonlinear heat equa- 
tions of the type 

d t u = Au + F(u, Vu, VVu). (1) 
and proved that, for a large class of nonlinearities F, the solution behaved, as t — ► oo, 



ed 



u(x,t)^At^r(^=) (2) 

with /* Gaussian, = e 4 . The asymptotics was thus proven to be "universal", 

i.e. independent on the initial data u(x, 0) and the nonlinearity F, within some class. 
The only dependence of them occurs in the Constant A. We explained this universality, 
following the work of Barenblatt [0] and of Goldenfeld et al @, §], in terms of the 
Renormalization Group, which also was an important ingredient in the proof (see Section 
3). 

The assumptions for the data in fìj was smallness in a suitable norm implying falloff 
at infinity (for a similar approach with initial data not decaying at infiriity, i.e. the 
formation of fronts and patterns, see || and for large data blowing up in a finite time, 
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see 0). The assumptions for F excluded terms sudi as u p , for 1 < p < 1 + -, where n 
is the dimension of space (p = 1 + - was considered and in this case the Constant A in 
(TJ) is replaced by A(\ogt)~%). 

Here we want to consider the asymptotics of ([TJ) in the presence of these "relevant" 
(in the Renormalization Group terminology) terms in F. Although again rather general 
F may be considered, we formulate below our results for the nonlinear heat equation 

d t u = Au- \u\ p -\ (3) 

where u = u(x,t), x E M n , 1 < p < 1 + - , and we discuss generalizations in Section 3. 
Equation has a one-parameter family of scale invariant solutions (see Section 2) 

u(x,t) = F"^f 1 {xt~^) (4) 
where 7 > j p > 0, provided 1 < p < 1+| and / 7 is non-Gaussian, having the asymptotics 

/ 7 (x) ~ \x\~^ (5) 
as | re | — ► 00 if 7 > 7 P , while, for 7 = j p , it decays at infinity as 



f lp (x) ~ |x|j>-i n e . (6) 

We prove here that these solutions are stable in two senses: first, there exists a ball 
in a Banach space of initial data such that the corresponding solutions tend, in the 
appropriate norm, as t — > 00, to (§]). Secondly, any initial data satisfying a suitable 
positivity condition will give rise to a solution again tending to (f|). 

More precisely, let q > and consider the Banach space B of L°° functions h 
equipped with the norm (with some abuse of notation!) 

11/11100 = esssup + (7) 

We consider the initial data (taken at time 1 for later convenience) 

u(x, 1) = f-y(x) + h(x) (8) 

with h G B. We prove the 

Theorem Let 1 < p < 1 + -. There exist e > 0, C < 00 and /i > suc/i t/iat, i/ the 
initial data u(x, 1) of (R|) g'wen (^) w//i h E B and satisfies either 

\\h\\oo < £ 

or 

h(x) > 

(a. e.) then, (^|) has a unique classical solution and, for ali t, 

||^«(^,t)-/ r (.)|| 0O <Ct-' 1 
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2 Proof 



Before going to the proof of the Theorem, we will briefly discuss the scale invariant 
solutions (f|). These are given by / 7 (x) = 7 (|x|) and 7 solves the ordinary differential 
equation 




- <P P = (9) 

\ '/ * j v — ± 

for j] = \x\ e [0, oo[. 

The theory of positive solutions of @ has been developped in || 0, • The main 



result is that, for any p > 1, there exists smooth, everywhere positive solutions, 7 , of 
@ with 7 (O) = and 7 (O) = 7 for 7 larger than a certain criticai value 7 P (but not 
too large). Actually, for p < 1 + -, j p > while 7 P = for p > 1 + -. The decay at 
infinity of these solutions is given in ([5], 0). 

The existence of a criticai 7 P can be understood intuitively by viewing (§) as Newton's 
equation for a particle of mass one, whose "position" as a function of "time" is 4>{rf). 
The potential is then U{4>) = 2 (p-i) ~ ^fT anc ^ the "friction terni" + |) 0' depends 
on the "time" rj. Hence, if 7 (O) = and 7 (O) = 7 is large enough, the time it takes to 
approach zero is long and, by then, the friction term has become sufficiently strong to 
prevent "overshooting" . However, as p increases, the potential becomes flatter and one 
therefore expects j p to decrease with p. 

Given the initial data (§), it is convenient to rewrite (|3]) in terms of the variables 
£ = xt~2 and r = logt; so, define v(£, r) by: 

u(x,t) = t~^{f^(xt~^) + v(xt-^ìogt)) (10) 

where now 

v(t,0) = h(£). (11) 
Then, (|3]) is equivalent to the equation 

d T v = Cv- (|/ 7 + vr\f y + v)- - pf^v) =£v + N(v) (12) 

where we used the fact that ([|) solves (Q) and gathered the linear terms in 

C = C + Vry, 

with 

£ = A + Ì -V+— (13) 
2 p — 1 

and 

To prove that the solution t~~ (f 7 (xt~^)) is stable means to find a class of initial 
data v(£, 0) such that the corresponding solution of (0) goes to zero as r — » 00, in a 
suitable norm. 
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The Theorem of Section 1 reads now in terms of v as 



Proposition 1. With the assumptions of the Theorem, with initial data (11) ,has 
a unique classical solution and 

IK-,T)|| 00 <Ce-'* T 



The main input in the proof is the following estimate on the semigroup e TC : 

Proposition 2. The operator e rC is a bounded operator in the Banach space B, and its 
norm satisfies 

\\e TC \\< Ce~^ T 

for some fi > 0, C < oo. 

There are two important ingredients in the proof of Proposition 2. The first is the 
fact that e rC is a contraction in a suitable Hilbert space of rapidly decreasing functions. 
To see this, note first that Co is conjugate to the Schròdinger operator 

e 2 e 2 f 2 n 1 

eV£ e-V = A-f--^ + ^- (15) 
lo 4 p — 1 

i.e. the harmonic oscillator. Thus Co is self-adjoint on its domain T>(Cq) C L 2 (R™, dfi), 
where 

d(J>{£) = e 4 d£. 

Cq has a pure point spectrum — f — y|m = 0, 1,...} and the largest eigenvalue 
-^t — f is positive ifl < p < 1 + -. Thus e rC ° is not contractive and, for e T ( Co+v ~<> 
to contract, we need to use the potential in a non-trivial way (this is the reason why 
l<n<l + -is harder than the p > 1 + - case) . 

Remarkably, it is possible to prove that C < —E < without a detailed study of the 
function / 7 , but only using equation @. We have the 

Lemma 1. The operator e rC is a bounded operator in the Hilbert space L 2 (R n , dfi) and 
its norm satisfies 



\e rC \\ < e~ Er 



for some E > 0. 



Proof. Since Vy is bounded, C is self-adjoint and, as for Co, its resolvent is compact 
and, therefore, its spectrum is pure point. Let — £7 7 be the largest eigenvalue. 

First note that — E 1 < ~E lp . Indeed, this holds since Vy < V lp1 because / 7 > f lp , 
which in turn follows from the fact that 7 P is the smallest allowed value of 7 (O) = 7 in 
(||), and that two solutions of (Q), both with initial conditions 4>L(0) = 0, will not cross. 
Hence it suffices to prove the claim for 7 = 7 p . Let us write E = E lp . 

Next we note that, by the Feynman-Kac formula [13], e rC has a strictly positive 
kernel; indeed, since — C < V y (£) < 0, we have 

e- rC e^°(e,0 < e r£ (^,n < e TC °(^0 (16) 
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and e TC ° is explicit, see flT8|) below. Thus, by the Perron-Frobenius theorem |TT] 



£ has a unique eigenvector fi with eigenvalue —E and Q can be chosen to be strictly 
positive. That —E is strictly negative, is now shown by the following argument. 
Assume for a moment that / = f Jp G T>(£), and write © as 

Cf = -( P -i)F. 

So, 

(fW) = -( P -i)(fì,f) 

where (•,•) denotes the scalar product in L 2 (R n ,dfi). By the self-adjointness of £ and 
the definition of Q, i.e. CQ = —EQ, we have 

- E =- (! '- 1) 7a7r <0 

since £7 and / are strictly positive. 

Finally, to see that f 7p G T>{C) we use the fact that, since f 7p decays at infmity like 
@, we have, for any 6 > 0, and any n G N, 

|V^/ 7p (0l<C(5,n)e-^^ 2 . (17) 

This follows easily from @ and the differential equation (||), and implies that f Jp G 
V(C )=V(C). ^ □ 

Notice that functions in L 2 (R, n ,dn) have essentially a Gaussian decay at infmity, 
which is much faster than what is allowed in our Banach space B, see (H). This brings 
us to the other cruciai ingredient in the proof of Proposition 2, which is that é rCo 
contracts functions in B pointwise for £ large. This follows from the explicit formula 
(Mehler's formula PI): 

(e rCo mO = (4tt(1 - e-))" V(iVt) exp {- ^^'^ ) ( 18 ) 
Hence, if a function v satisfies 

wo\<c(i+\m~\ (19) 

for some Constant C, we have 

Ke^MlOI^C'ei^a + lel^)- 1 (20) 

for |£| large enough (of order a/t) and another Constant C Hence, the operator é rC ° 
contracts, for large |£| and large r, any function that decays as in © with g > By 
(p~6|) , we see that C behaves similarly. 

The idea of the proof of Proposition 2 is the following. For |£| small (EO) seems to 
expand by e~: the potential y is important in this region and we want to use the 
information we obtained in the Hilbert space, Lemma 1 (recali that these functions have 
rapid decay, so that this bound should be used to capture the contraction only in the 
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small £ region). For large £ , we shall use (20). This small-large £ interplay is however 
slightly subtle, and we need to resort to an inductive argument to control the large r 
behaviour in Proposition 2 (this is actually just the Renormalization Group idea applied 
to the linear problem). 

Proof of Proposition 2. It is convenient to introduce the characteristic functions 

Xs = < p) 

Xe = > p) 

where p will be chosen suitably below. The properties of C that we need are summarized 
in the following 

Lemma 2. There exist constants C < oc, E > 0, and 5 > 0, such that 
i) For g E B, 

Ile^Hoo^Cei^lMU. (21) 

ii) For g G L 2 (R n , dpi), 

||e £ <?||oo<C7|kl|2, (22) 
where || • H2 is the norm in L 2 (R n , dfi) . 

iii) For g such that Xs9 £ L 2 (R n , dfj), 

c -i- 
\\Xee p Xs9\\oc < e 5 \\Xs9\\2, (23) 

for p large enough. 

iv) For g G B, 

HX^IU < e-^I^IU, (24) 

for p large enough. 



Let II (7 II oo = 1. Given Lemma 2, we set r n = np, and prove inductively that there 
exists a > such that v{r n ) = e TnC g satisfies, for p large, 



\Xsv(r n )\\ 2 + \\xsv(t, 



< e~e~ an 

ni II 00 _ ; 



(25) 



and 

||x^(Tn)||oo<e- an . (26) 

Proposition 2 follows from (^,^), by taking p = - (for times not of the form r = np, 
use (pi])). The bounds ( p5] , p6] ) hold for n — 0, for p large enough, using ||g||oo = 1 and 
the obvious inequality 



\Xs9h < P 2 e » 



(27) 
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So, let us assume (|25| , p6|) for some n > and prove it for n + 1. Let t> = v{r n ) and 
write 

« = XsV + X£V = v s + v e . 
We have from Lemma 1 and (25) 

„2 



\e pC v s \\ 2 < e- pE e^e~ an 



and, combining ( p^) and Lemma 1 



|e p£ ^||oo < C||e (p ~ 1)£ iU 2 < Ce~ (p - 1)£; e^e- an . 



Finally, from ( |2Ì"D and ([26|), we have 

||e p£ ^||oo < Ce^e" 

and, from this and (^7]), we get 
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(29) 



(30) 



(31) 



Combining (p8|-|3ì"|) , one gets fl25|), with n replaced by n + 1 for p large enough and 
a small. On the other hand, (|26|) , with n replaced by n + 1, follows immediately from 
PIPI) and (PSIEI), taking a < Sp. □ 



We are left with the 

Proof of Lemma 2. Part (i) follows immediately from flTE| ) and ([ZO]). 

For (ii), we use Schwartz' inequality and the bound 



sup j \e c (£,e)\ 2 d£ <oo, 



(32) 



which follows from ^ |X6|) and (|T8|) . 

For (iii) proceed as in (ii) by using Schwartz' inequality, but replace ( |3"2"D by 



sup 

\Z\>p 



e pC (c,o\ 2 x(\e\<p)de) <e-^ 



(33) 



which again follows from (|TBp and flTBp (we can replace | in ( |3"3"| ) by | — 5 for any 5 > 0, 
if p is large enough). 

Finally, (iv) follows from (^) and q > □ 

Proof of Proposition 1. The proof is straightforward, given Proposition 2. We 
consider the integrai equation corresponding to fll2|): 



(34) 



( T ) = e r£ /i+ / dse (r - s)£ iV(v(s)) = 5(v,r) 



with ufrì = v(-, t). 
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First, let H^lloo < e; (@) is solved by the contraction mapping principle in the Banach 
space B of functions v = v(^, r), where v(-, r) G B for r G [0, oo), and where the norm 
is 

||u|| = sup \\v(-,T)\\ 00 e Tfl 

rG[0,oo) 



We show that S defined by fl34|) maps the ball £>o = {v G B| ||i>|| < Ce} into itself, 
for a suitable Constant C, and is a contraction there . This follows, since we get, using 

\N(v)\ < C'\v\ p 
where p = min(2,p) > 1, and C is a Constant; therefore 

\\N{v{s))\\ 00 <C>\\v{s)\\l<C>é>e-^ (35) 

and so, Proposition 2 gives the claim for e small. The proof that S is a contraction is 
similar. Finally, from the Feynman-Kac formula and the smoothness and boundedness 
of V one deduces that e r£ (£,£ / ) = e T£o (£, £,')K(£, £') where K is smooth and bounded. 
The regularity of the kernel of é rC (for short times, it behaves like the heat kernel), 
implies that the solution of (34) is actually the unique classical solution of ( |T2"| ) (for 
details of such arguments, see 0). 

Let now h > 0. It is standard that equation ([3]) has unique solution for ali times, 
with positive initial data such as ours |J. So, equation (|i~2]) also has a unique solution 



for ali times. Moreover, by a comparison inequality (v — is a solution of (|ì~2|)), we have 
v(£, T ) > for ali times. Thus, in (|34]) , the second term is negative (flì~2|) shows that 
N(v) < 0, and (|T6| , 18) show that the kernel of e rC is positive). Therefore, we have the 
pointwise inequality 

v(t,r)<(e TC v(0))(0 

and Proposition 2 proves the claim. □ 



3 Extensions and concluding remarks 



1. In Theorems 1 and 2, we could use a norm defined by 

\\g\loo = esssup |#(£M£)|. 



(36) 



where w is any bounded positive function decaying at infinity faster than |£| - ~. Indeed, 
the norm (7) was used only in the proof of (20). However, we would not necessarily get 
exponential decay of v as a function of r. 

2. In |l0l and also in JT^j, results similar to ours were obtained on the stability of 
the self-similar solutions for p < 1 + ~ and 7 = 7 P . However, our results apply to a 



ball in a Banach space (and to any 7) while in [10| the initial data is assumed to satisfy 
a pointwise inequality (but not to be small). This is similar to the h > case in our 
theorem, but with another inequality. For the first part of our Theorem, u(x, 1), given 
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by (8), does not even have to be positive. On the other hand, very general results ori 
the stability of the self-similar solutions, decaying as in (5), were obtained in ||13|| , for 



p > 1 + -. Basically, it is shown there that any positive initial data decaying at infinity 

2 

like |£| p- 1 will give rise to a solution whose asymptotic behaviour in time is given by 

3. With little extra work, the small ||u|| part of the Theorem generalizes to more 
general nonlinearities, e.g. equations of the forni 

d t u = Au - \u\ p ' l u + F(u, Vu) (37) 

whereby we need to add to (|12|) the term 

F T (v,Vv) = e^F (e~^(f 7 + v),e~^^ T V(f 7 + v)\ . (38) 

In order to define "small" initial data, we introduce the Banach space of C 1 functions 
with the norm: 

||fc|| = ||/i||oo+ max \\dih\\co (39) 

l<i<n 

where ||/i||oo is defined in (7). 

We assume that F in ( |37|) is C 1 and satisfies: 

\F(a,b) \ + \ad a F(a,b)] + max |6id 6 ,F(a, 6)1 < X\a\ Ql ( max |ò 4 |Y 2 (40) 

l<i<n 1 / 

for \a\, \bi\ < 1, where 

v+ 1 

qi + ^-g 2 > p (41) 

and A is taken small. Using (40, 41), and the fact that / 7 , dif 7 , belong to the space B, 
one gets 

Valico <CAe- 5s 

for some 5 > 0, which is like the RHS of (35) for A small. The last two terms in the 
LHS of (40) are used to prove that the F term defines a contraction. The only extra 
difficulty is to show that the solution of the fixed point problem (corresponding to (34)) 
is a classical solution. However, using the regularity of the kernel of e TC discussed in 
the proof of Proposition 1, one shows first that v is C 2 ~ a for any a > (its derivative 
is Hòlder continuous of exponent 1 — a); then, one uses that information, the integrai 
equation and the regularity of the kernel to show that v is actually C 2 . See [0 for details. 

4. Finally we want to end with some comments on the RG picture behind these 
results. In []l] the RG map TZl, for L > 1 was defined in a suitable Banach space 
of initial data f(x) = u(x, 1) and a suitable space of nonlinearities F (taken to be 
holomorphic functions in U). TZl consisted simply of solving ([!]) up to the finite time 
L 2 and performing a scale transformation on the solution and on F: 

n L (f,F) = (f L ,F L ) (42) 
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where f L (x) = L n u(Lx,L 2 ) and F L (u,v,w) = L 2+n F(L~ n u, L~ n ~ x v, L~ n ~ 2 w). This 
scaling assures the semigroup property lZ L k = lZ k L on a common domain. The limit (0) 
was then shown to follow from 

n h L (f,F)^(Af*,0) (43) 

as k — > oo, where (Af*,0) is a one parameter family of Gaussian (f*(x) = e~~) fixed 
points for TZl. Universality, i.e. independence on / and F was then explained in terms 
of a dynamical systems picture: if (/, F) lie on the stable manifold of the line of fixed 
points, ali the corresponding equations and data have the same asymptotics. 

The Theorem of Section 1 is a statement on the stability of a family of non- Gaussian 
fixed points of the RG (which are, moreover, non perturbative, i.e. we do not use any 
"e-expansion" ) . Equation (1) is invariant under the scaling u — > ul with 

2 „ 

UL{X,t) = Lp~ 1 u(Lx, L t) 

which suggests setting 

f L (x) = L^u{Lx,L 2 ) 

and a correspondig definition of F L . Then TZ L {f 1 , F*) = (f 7 ,F*) for F*(u) = —u p , and 
our Theorem constructs the stable manifold of this fixed point (in the / variable; for 
F, see the previous remark). The reason one needs an iterative approach to the limit 
(|43|), i.e. one controls the iteration of TZl rather than 1Z\ directly, is the existence of 
the neutral direction: A is a nontrivial function of the data and of the equation. Here 
there is no neutral direction (this is the content of Lemma 1), and no iteration is needed 
(although we needed to resort to one in the proof of Proposition 2! But that was for a 
different purpose, namely the analysis of the linear operator é rC ). 

Combining our results with those of [Q , we obtain the following p-dependence of the 
asymptotics of the solution of 

d t u = d 2 u — u p , 
for a suitable class of initial data. For p > 1 + -, 

r ri' 

ui-t^^^At-^n) 

where /*(£) = e - ^ 2 / 4 is independent of p; the prefactor A depends on p and on the initial 
data. We have 

J |w(a;,t)|da; = 0(1) 

as t — > oo. For p — 1 + -, 

i.e. it is as before, but with a logarithmic correction and 

f \u{x,t)\dx = 0{{logt)- 1/2 ). 
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For Kp< 1 + f, 

u(-& 2 ,t)~t-èiM-) 

where / 7 is non-Gaussian, varies with p but is independent of the initial data (satisfying 
the hypotheses of the Theorem). Moreover, 

J \u(x,t)\dx = O -> 

as t — > oo. 
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